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We study the properties of the free energy of infinitely heavy quark anti-quark pair in SU(2) gauge 
theory. By means of lattice Monte Carlo simulations we calculated the free energies in the singlet, 
triplet and color averaged channels, both in the confinement and in the deconfinement phase. The 
singlet and triplet free energies are defined in Coulomb gauge which is equivalent to their gauge 
invariant definitions recently introduced by Philipsen. We analyzed the short and the long distance 
behavior, making comparisons with the zero temperature case. The temperature dependence of the 
£f>^ • electric screening mass is carefully investigated. The order of the deconfining transition is manifest 

' in the results near T c and it allows a reliable test of a recently proposed method to renormalize the 

, Polyakov loop. 
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I. INTRODUCTION 



The study of the free energy of a static heavy quark-antiquark pair in a medium of color charges at some temperature 
T has recently become a hot topic of statistical QCD PJ-3- The main reason is the fact that such free energy is 
related to the potential between quarks and antiquarks, which is of fundamental importance both for understanding 
deconfinement and in heavy quark phenomenology at finite temperature 0-0 The study of static quark free energies 
(Polyakov loop correlators) is also important for constructing effective theories at the deconfinement transition |7J 

The presence of the medium affects the quark-antiquark interaction in a non trivial way. Such modifications of 
interactions is usually studied in terms of the free energy of a static quark-antiquark pair separated by some distance R. 
So far most studies concentrated on the behavior of the static quark-antiquark free energy at large distances, R 3> 1/T 
HI ; below the deconfinement temperature T c such behavior is characterized by a temperature dependent string 
tension, above T c by exponential screening which is governed by a temperature dependent (Debye) screening mass. 
qj*' However, it turned out that the free ener gies of static quark-antiquark exhibit a quite complex temperature behavior 
already at short distances R < 1 /T 0, llfl fTTL Il2| . Furthermore, the detailed study of the free energy at short 
J* 1 distances allows to define a renormalized order parameter [3]. As far as the physics of heavy quarkonia is concerned, 
. ^ \ the detailed structure of the static quark-antiquark free energy at short distances is even more important than its 
large distance behavior. Though the large distance behavior of the color averaged potential was extensively studied in 
Refs. 0,111! the problem was not completely settled. This is partly due to finite size effects and very large statistics 
needed in such studies. 

In this paper we want to study various features of the heavy quark free energy in SU (2) gauge theory at finite 
temperature. The study of SU(2) lattice gauge theory has at least two advantages: the simulations are not so time 
consuming as in ST/(3) and the deconfinement transition is second order, which has interesting consequences on the 
behavior of the free energy near the critical temperature T c . In general, the heavy quark free energy depends on the 
color channel one considers; for a complete analysis we investigated the singlet, the triplet and the color averaged 
channels. In particular the study of the color singlet free energy is of interest because it is the most relevant quantity 
as far as the physics of heavy quarkonia at finite temperature is concerned. We examined both the short and the 
long-distance behavior of the free energies, below and above the deconfinement temperature. We devoted a special 
attention to the issue of screening, in that we determined the Debye masses at various temperatures. 

The rest of the paper is organized as follows. In section II we define the free energy of a static quark-antiquark 
pair in color singlet, color triplet and color averaged channels and discuss the choice of the simulation parameters. 
The basic features of the static quark-antiquark free energies are also discussed there. In section III we present 
the numerical results below T c . Section IV deals with free energies above deconfinement and determination of the 
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screening masses. In section V we define the renormalized Polyakov loop for SU(2) gauge theory following Ref. 0]. 
Finally section VI contains our conclusions. 

II. FREE ENERGIES IN SU(2) GAUGE THEORY 

On the lattice the free energy of a static quark-antiquark pair in the gluonic medium is determined by correlation 
functions of temporal Wilson lines L, 

N T -1 

L(R)= ft U (R,t), (1) 

T = 

TrL is also referred to as the Polyakov loop. Following Refs. ^| 0] we introduce the color singlet and triplet free 
energy of a static quark-antiquark pair 

e -F l{R ,T)/T + c = }.(Tr(L(R)tf(0))), (2) 

e -F 3 (R,T)/T+c = l(TrL(R)TrL^{6)) - -(TrL(R)LUo)) (3) 
3 6 

as well as the color averaged free energy defined by 

e -F^ { R,T)/T+c = l(TrL(R)Trtf(0)). (4) 
The latter can be written as a thermal average of the free energies in singlet and triplet channels, hence the name, 

e -F avg (R,T)/T = \ ( ,-F 1 (R.T)/T + ^_ e -F 3 (R,T)/T _ ^ 

The normalization constant C can be defined in different ways. In the deconfmed phase it is customary to set 
C = In \ (^TrL}\ . Another possibility is to fix it by normalizing the singlet free energy to the zero temperature heavy 
quark potential [3j . 

The main problem with the definitions of the singlet and triplet free energies ©-© is that these definitions are not 
gauge invariant, as the Wilson line is not a gauge invariant quantity. The only manifestly gauge invariant quantity is 
the color averaged free energy. This is the reason why singlet and triplet free energies were not studied in much detail 
so far. It was recently showed by Philipsen that gauge invariant definitions of the singlet and triplet free energies can 
be achieved by replacing the Wilson line in Eqs. 101, © by a gauge invariant Wilson line defined by 

L(R) = n\R)L(R)n(R) (6) 

The SU (2) matrix f2(i?) is constructed from eigenvectors of the spatial covariant Laplacian (see Ref. 0] for further 
details). Furthermore it was shown that this definition is equivalent to the definitions of the singlet and triplet free 
energies in Coulomb gauge. Since the determination of eigenvectors of the covariant Laplacian is computationally 
very expensive we fix the Coulomb gauge to calculate the singlet and triplet free energies. 

In our numerical investigations we use the standard Wilson action. In order to get control over finite size effects, 
which become important in the vicinity of T c , we have performed simulations at several different volumes. As we also 
want to investigate the short distance behavior of the free energies simulations were performed for N T = 4, 6, 8. To fix 
the temperature scale we have used the non-perturbative beta function of Ref. [l5|]. We will also use T c /y/a = 0.69 
|l6j |. with a being the zero temperature string tension. The lattice volumes and the gauge coupling along with the 
corresponding temperatures used in our simulations are summarized in Table 1. 

Calculations of the zero temperature potentials and of the free energy of a static quark-antiquark pair show violation 
of rotational symmetry at short distances. Since we are also interested in the behavior of the potential at short 
distances, we should try to remove these lattice artifacts. Following Ref. 0] we replace Fi(r) by Fi(ri) where 
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2.2900 0.97 32,60 

2.2930 0.98 60 

2.2960 0.99 32,60 

2.2975 0.995 32 



2.3019 1.01 32,48,60 

2.3077 1.03 32 

2.3134 1.05 16,32,48 

2.3272 1.10 32 

2.3533 1.20 16,32 

2.3776 1.30 32 

2.4215 1.50 16,32 

2.5118 2.00 16,32 

2.6431 3.00 32 

2.8800 6.062 32 

2.9766 8.000 32 

3.0230 9.143 32 

3.2190 15.87 16 



TABLE I: Lattice volumes and gau ge couplings /3 = 4/g 2 used in our simulations. The values of T/T c were obtained using the 
non-perturbative beta function flql . 



In this way we replace the lattice separation by the separation rj which corrects for the tree level artifacts in the 
Coulomb potential calculated on the lattice. When presenting the data on the free energy we will always do this 
replacement unless stated otherwise. 

Let us now present some general features of our findings. First we have performed simulations at a fixed lattice 
spacing, corresponding to the gauge coupling (3 = 2.5 at iV r = 12 and N T — 6 corresponding to temperatures T — 0.6T C 
and T — 1.3T C , respectively. The results are shown in Fig. Q](left). At this value of the gauge coupling the ground 
state static quark-antiquark potential as well as the first two excited potentials at T = have been calculated [18j and 
we show them in Fig. ^together with our data. The singlet free energies at short distances do not differ from the zero 
temperature potential and temperature dependence shows up only at r^/a > 0.5. On the other hand we find that the 
triplet free energy is considerably smaller than the first excited potential at T — 0. At small distances it was shown 
that the excited potential coincides with the perturbative octet (triplet) potential up to a non-perturbative constant 
|l9j | . No such constant is expected in our definition of the triplet free energy which we expect to coincide with the 
perturbative one at short distances. At very short distances also the triplet free energy is temperature independent. 

In most of our calculations we have varied the temperature T by varying the lattice spacing a (i.e. the gauge 
coupling (3 = 4/g 2 ) for fixed temporal extent N T . As both the free energy and the T = potential contain a lattice 
spacing dependent additive constant C (c.f. Eqs. 1215(1 . some normalization prescription should be introduced in order 
to compare the free energies calculated at different temperatures. To do so we assume the following form for the zero 
temperature potential 

rr , w ,- 0.238 r- 0.0031 

V{r)yfd= ^+rV^4-—— - (8) 

This form was obtained in Ref. [2(3 by fitting the lattice data on the T = potential for r^fa > 0.063 at /3 = 2.85 
apart from the constant which we have omitted. The violation of rotational invariance on the lattice was taken into 
account in the fit procedure. Thus equation JSJ defines our convention for the continuum zero temperature potential. 
In some cases we need the T — potential at distances r^/a < 0.063. In this case we use the 3-loop perturbative 
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FIG. 1: (Left) Free energies of a static quark-antiquark pair as a function of the distance r^/a at T — 0.6T C (open symbols) and 
T — I.3T C (filled symbols) corresponding to coupling /3 = 2.5 Also shown there is the T = potential and its first excitation 
(open and filled lower triangles correspondingly connected by lines) at the same (3 value. (Right) The free energies are here 
calculated at N T = 4; the temperature values are T = 0.9T C and T = I .3T C , respectively. The triangles, squares and circles 
indicate the singlet, the triplet, and the averaged free energies, respectively. 



potential calculated in qq scheme [2lJ normalized to smoothly match the form defined by Eq. JSJl at r^fa = 0.07. We 
have also checked that the difference between the 3-loop and 2-loop results is negligible for our purposes. In what 
follows the normalization constant C will be chosen such (unless stated otherwise) that the singlet free energy matches 
the T = potential at the shortest distance (r/a = 1). In Fig. ^ (right) we show our results on static free energies 
calculated for N T — 4 and using this normalization convention. The free energies in the deconfined phase reach the 
same value at large distances (see Fig. P). This is to be expected as at very large distances due to screening the free 
energy of color charges should be independent of their relative color orientation. What is more interesting is that at 
large distances the gap between the triplet and singlet free energy vanishes as well in the confinement phase though 
it is non-zero for the T — case. This issue will be discussed more in detail in the next section. 



III. RESULTS IN THE CONFINEMENT PHASE 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 



FIG. 2: The color singlet (left) and triplet (right) free energies at various temperatures in the confinement phase calculated for 
N T = 4. The normalization constant C was chosen such that the singlet free energy matches the zero temperature potential 
(solid line) at the shortest distance (see text). 



In this section we are going to present our numerical results in the confinement phase. In Fig. [3we show the singlet 
and triplet free energies with the normalization described in the previous section. One can see that the temperature 
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FIG. 3: (Left) Color averaged free energies at various temperatures in the confinement phase, (right) String tension as a 
function of T/T c , crosses correspond to 32 3 x 4 lattice, open squares to 60 3 x 4 lattice. 



dependence of the singlet free energy is only visible at distances r^fa > 2. The triplet free energy is also temperature 
independent at small distances; thermal effects become visible at r^fu > 1.4. We also notice that there is a slight 
enhancement of the singlet free energy over the T = potential in the interval 1 < r-s/a < 2. A similar enhancement 
was observed also in the case of 2 + 1 dimensional SU(2) gauge theory at finite temperature [2j as well as in preliminary 
SU(3) calculations 4]. As the color averaged free energy is a thermal average of singlet and triplet free energies, it 
would have a non-trivial temperature dependence even if the latter were temperature independent. This temperature 
dependence is larger the smaller the gap between the triplet and singlet contributions. In general one can say that the 
temperature dependence of the color averaged free energy at short and intermediate distances is mostly determined 
by the value and the temperature dependence of the color triplet contribution. If this continues to be true in full 
QCD, then some of the conclusions of Ref. (where the color averaged free energy was related to the meson masses) 
should be revised. 

At large distances the color singlet, triplet and averaged free energy reach a common value which can be parametrized 
by a form 

Ft{r, T)|rT>i = cr(T)r + A(T) InrT + B(T) (9) 

In Fig. |3| we show the color averaged free energy as well as the string tension obtained from it using the fit to Eq. 
[5] Because the deconfinement transition is of second order, the inverse correlation length (string tension) vanishes 
at T c . As a result we have large finite size effects close to T c . Indeed we found that on the 32 3 x4 lattice the string 
tension vanishes around 0.97T C . For this reason we performed simulations on the larger 60 3 x4 lattice close to T c . We 
have seen that the values of a(T) we have obtained from our calculations on 32 3 x 4 lattice agree with the results 
from the 60 3 x4 lattice up to 0.92T C . Above such temperature we adopt the results from the larger lattice. The final 
situation is illustrated in Fig. El (right), where the squares indicate the results from 60 3 x4, the crosses the ones from 
32 3 x4. It is well known that tr(T) vanishes near T c according to the power law cr(T) oc (T c — T) v ', where v = 0.63 
is the 3D Ising exponent for the correlation length. Therefore we tried to fit our data point by using the ansatz 
ct(T) = a(T c — T) u [\ + b(T c — T) 1 ' 2 ], with v — 0.63. The best fit curve is shown in our plot and it reproduces very well 
our data. The values of cr(T) found by us are considerably smaller than those obtained in Ref. This is probably 
due to the fact that the lattice volumes used in Ref. || were considerably smaller than ours. 



IV. RESULTS IN THE DECONFINEMENT PHASE: SCREENING 



We start our discussion of the numerical results in the deconfined phase with Fig. 01 where we show the singlet free 
energy in units of ^fa normalized to the zero temperature potential at the shortest distance. As one can see the singlet 
free energy saturates at large distances, while at short enough distances it is temperature independent and coincides 
with the zero temperature potential. The distance in physical units, at which the temperature dependence enters, of 
course strongly depends on the value of the temperature, the higher the temperature the shorter is the distance where 
effects of the medium become visible. Also the distance where the singlet free energy saturates strongly depends on 
the temperature, it is getting larger as we approach T c . Close to T c screening enters only at distances r^/a > 1. This 
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feature of the free energy is reflected in the temperature dependence of the screening masses which will be discussed 
below. Another interesting feature of the singlet free energy is that the value of the plateau of the free energy decreases 
with the temperature. Such a behavior of the singlet free energy was observed for SU(3) gauge theory in Ref. 0], 
where it was also argued that the reason for this is the presence of the entropy contribution. 

For the further discussion of the results in the deconfined ph ase, especially for making comparisons with perturbation 
theory, it is more convenient and in fact customary (cf. |lL to choose the renormalization constant C in Eqs. 
@, © to be C = In \(^TrL)\ 2 . Furthermore the large distance behavior of the free energies should be discussed 
separately from their short distance behavior, where one would expect perturbation theory to work. In general it is 
expected that perturbation theory breaks down at distances r > l/g 2 T [23, with g being the gauge coupling constant. 
As for the physically interesting temperature range one always has g ~ 1, perturbation theory may be applicable at 
distances rT < 1. One of the predictions of perturbation theory is that — 3F^(r, T)/Fi(r, T) ~ 1, for any r. In Fig. 
[3] we show this ratio for different temperatures. This ratio appears indeed to be constant for T > 1.5T C but always 
smaller than 1 , even for the highest temperature we considered (16T C ). Similar results were found in Landau gauge 
inRefs. ESEl. 

High temperature perturbation theory predicts that the color averaged free energy has the form [l3l |2^ | 



F avg (r,T) _ 3 ff 4 2m D0 r / in \ 

T 32(4^rT) 2 1 ' 

at distances r > 1/T, with m D0 being the leading order Debye mass mno = \/2/3gT. At distances r < 1/T the 
simple form 1)11) [) is no longer valid, though the l/r 2 -like behavior is still expected due to cancellation between singlet 
and triplet contributions |l3j . Therefore we define the so-called screening function S(r,T) by the formula 



F avg (r,T) _ 3 1 



; S(r,T) (11) 



T 32 (rT) 

In Fig. 0we show the numerical results for the square root screening function. As one can see from the figure, at high 
temperatures (T > 2T C ) the screening function shows a mild r-dependence, which implies that the color averaged free 
energy behaves like 1/r 2 . The screening function S(r, T) decreases with with increasing thenperature which one would 
expect if S(r, T) ~ g 4 (T). At temperatures closer to T c the screening function decreases at small distances. Obviously, 
this behavior has nothing to do with screening and signals the breakdown of the high temperature expansion. As we 
approach T c Fi^/T is no longer small at small distances as shown in Fig. 6 and therefore the exponentials in Eq. (J5J 
cannot be expanded. As a result of this, the cancellation between the singlet and triplet free energies no longer holds; 
moreover, since the singlet free energy is negative and the triplet one is positive (when working with normalization 
convention C — In |( ^TrL) | 2 ) the color averaged free energy is dominated by the singlet contribution and behaves like 
1/r at small distances. We also note that the free energies calculated for different N T (different lattice spacing) at 
3T C agree reasonably well. Similar agreement between the results calculated for different values of N T was observed 
for other temperatures too. 
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FIG. 4: The color singlet free energy at various temperatures in the deconfinement phase calculated for N T = 4 (left) and 
N T = 8 (right). The solid line is the T = potential. 
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FIG. 5: Free energy triplet/singlet ratio at various temperatures. 
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FIG. 6: The square root of the screening function (left) and the singlet free energy in units of T (right) at different temperatures. 
The open symbols refer to the results for N T = 4, the full symbols to the results for N T = 8. 



Let us now discuss the large distance behavior of the free energies and the determination of the screening masses. 
We will restrict ourselves to the discussion of the color singlet and color averaged free energy as the color triplet free 
energy becomes very noisy at large distances and the present statistics does not allow to study it in detail. 

Contrary to earlier studies 0, 0, EI j where the screening masses were obtained using uncorrelated fit and the short 
and large distance behaviors of the free energy were not separate, here we use the correlated fit procedure of Ref. 
[25j . From Fig. 6 it is clear that in the region rT < 1 the color averaged free energy can be well described by almost 
unscreened l/r 2 -like behavior. Therefore this region should not be considered for the determination of the screening 
masses. In our procedure the fit interval was chosen so that the fit yields a reasonable x 2 / .o.F). Thus, unlike in 
the uncorrelated fit used in Ref. pj, there is no dependence of the screening masses on the fit interval. We determine 
the screening mass in the color singlet channel by fitting the data with a screened Coulomb (Yukawa-like) ansatz. In 
leading order perturbation theory, in fact, the most important contribution to the singlet free energy is given by the 
exchange of a single gluon. In the case of the color average free energy we used a more general fit ansatz |32j : 

*jp = -Aexp(- M i?) + S. (12) 

As possible values for the exponent d we took d = 1,2. The resulting values of the screening masses that we extracted 
are presented in Tables [n] and IIIII for the singlet and the averaged channels, respectively. For determination of the 
screening masses we mostly use 32 3 x 4 lattice. 

We found that for the color averaged free energy the fits with d = 1 and d = 2 are both good for all temperatures, 
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except near T c , where we got a reasonable \ 2 1 O^D.o.F) only for d = 1. From our fit analysis it then comes out that 
we cannot choose between the two cases, which gives a systematic error of about 30% on the screening masses. In 
order to eliminate this ambiguity, we have also calculated the plane-plane correlator, given by the formula: 

C PL (x 3 ) = (Tr L(x 3 )Tr L*(0)) - \{TrL)\ 2 (13) 

where L(x 3 ) = J2 Xl x 2 L{xi,X2,x 3 ). If the color averaged free energy has the form flT^} with d = 1, then Cpl{x 3 ) 
should fall off with the distance as a simple exponential, which allows a direct determination of the screening mass. 
The results we found are reported for comparison in Table ITTT1 We see that the values of the masses extracted from 
the plane-plane correlator agree in each case with the masses obtained from the point-point correlator when d = 1. 
Furthermore, we have analyzed the effective masses extracted from Cpl(x 3 ). They reach a plateau already at x 3 T ~ 1, 
which makes the presence of power-like prefactors in Cpl(x 3 ) at large distances very unlikely and thus implying that 
very likely d ~ 1. At high temperature dimensional reduction arguments suggest that d — 1, as the large distance 
behavior of any static correlators is governed by exchange of a bound state of the effective three dimensional theory 

For T — 1.01T C the results refer to the 48 3 x4 lattice instead of 32 3 x4. In this case, in fact, the mass changes 
appreciably when one goes to the larger 48 3 x4 lattice. This is probably due to the large correlation length close to 
T c . Further simulations on a 60 3 x4 lead to the same value of the mass we found on the 48 3 x4, which is then reliable 
as infinite volume limit at this temperature. 
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FIG. 7: Screening masses in units of temperature vs. T/T c extracted from the color singlet free energy (left) and from the 
color averaged free energy (right). In the left plot we have also shown the screening mass extracted from the static electric 
propagator [%j |26| . In the right plot we give the lowest Af screening mass from |28| as well as the screening masses of 
plane-plane correlators of Polyakov loops from [gsj (open triangles). 



In Fig. \7\ we show the screening masses extracted from the singlet free energy as a function of the temperature. 
There we also show the values of the screening masses obtained from the electric gluon propagator in Landau gauge 
[22Ll26l | at T > 1.2T C . As one can see they are compatible with the singlet masses we have found. This is to be expected 
in perturbation theory. But as g ~ 1 in the temperature interval considered and the screening mass governs the large 
distance behavior of the free energy, this agreement is quite non-trivial. In Fig. 0the color averaged screening masses 
are shown as well; we compared them with the lowest scalar screening mass (spatial glueball mass) obtained in 
Ref. . Dimensional reduction arguments [27j suggest that these masses should agree at high temperature and the 
figure seems to indicate that this is indeed the case. Very recently the plane-plane correlators of Polyakov loops were 
studied in Ref. [2^| for T = 1.1105T C and T — 1.227T C (we use the non-perturbative beta function to convert 
the gauge coupling of Ref. [2^ to temperature). In Fig. 0we show as well the corresponding screening masses (open 
triangles). Moreover we note that the color averaged screening masses obtained by us using fits with d = 2 agree with 
findings of Refs. H HI where the same exponent was used. The color averaged screening mass should go to zero 
when T^T C because it is just the inverse of the Polyakov loop correlation length, which diverges at T c . On the other 
hand there is no a priori argument for which the singlet mass should vanish at the threshold. Fig. [7| indicates that 
both masses become very small near T c . 
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V. THE RENORMALIZED POLYAKOV LOOP 



The last issue we would like to address is the renormalization of the Polyakov loop. It is known that if one takes the 
continuum limit at fixed temperature, the expectation value of the Polyakov loop vanishes, so that the usual definition 
does not really provide a physical order parameter for deconfinement. As we have already said at the beginning, the 
free energies on the lattice are always defined up to some renormalization constant. According to a recent work Q, a 
suitable choice of such renormalization constant can lead to a new definition of the Polyakov loop. If the constant is 
chosen such that the singlet free energy matches the zero temperature heavy quark potential at short distances, one 
can define a " renormalized" Polyakov loop L ren through the formula: 

L ren = exp(-F 00 (T)/2T), (14) 

where i ? o(T) is the asymptotic value of the singlet free energy at the temperature T (in fact, it does not depend on 
relative color orientation of the quark-antiquark pair, see above). We remind that we have renormalized the singlet 
free energy exactly in this way, so that, in our case, F^T) is nothing but the height of the plateau of the curves 
in Fig. 4. Practically we took F^T) = Fi(N a /2,T), which is the value of the free energy at the largest distance 
allowed on the lattice. In Fig. [SJwe plot the renormalized Polyakov loop as function of the temperatures. When 
normalizing the singlet free energy to the T — potential at the shortest distance {r/a = 1) we implicitly assume that 
there is no temperature dependence at this distance. This may not always be the case. Therefore we also calculate 
Foo{T) and the corresponding L ren by normalizing the singlet free energy to the T = potential as well at r/a = \/2. 
The difference in FodT) ( L ren ) arising from these two normalizations give us an estimate of possible systematic 
errors. When quoting the error on the renormalized Polyakov loop we always add quadratically the systematic and 
the statistical errors. At variance with SU(3), we now have a second order phase transition and a well defined scaling 
behavior of L ren at criticality. In order to check the scaling we have fitted the data on L ren using the standard ansatz 
in the interval T c < T < 1.5T C 



L ren (T) = c{T - T c f[l + b(T - T c )% (15) 

with the exponents (3 and uj fixed to their SU(2) values (3 = 0.3265 and uj = 1. The fit curve is shown in Fig. |H| 
(dashed line) and reproduces quite well the pattern of the data points. Considering only temperatures T < 1.1T C it 
is possible to fit the data with b = and approximately same value of c. 
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FIG. 8: Renormalized Polyakov loop as a function of the temperature in the entire temperature interval (left) and for T < 1.5T C 
together with the fit (right). The error bars indicate the combined statistical and systematic errors. 



VI. CONCLUSIONS 



In conclusion we have studied the quark-antiquark free energies in SU(2) gauge theory below and above the 
deconfinement temperature. We have found that the temperature dependence of the singlet free energy is much 
weaker than the temperature dependence of the color averaged free energy. Most of the temperature dependence of 
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the color averaged free energy is due to the presence of the color triplet contribution and its temperature dependence. 
If this will hold for real QCD it may have important consequences in the heavy quarkonia phenomenology at finite 
temperature. At large distances the free energy in color singlet and triplet channel converge to a common value. 

Above T c the color singlet free energy can be understood in terms of propagation of a non-perturbatively screened 
gluon. Through the detailed analysis of the point-point and plane-plane correlators as well as comparisons with other 
determinations of the screening masses we have established that the color averaged free energy is described by Yukawa 
law at large distances whereas at shorter distances it exhibits a more complex behavior. At low temperatures and 
short distances the color averaged free energy is dominated by the singlet contribution, while at higher temperatures 
it has l/r 2 T behavior and its temperature dependence is qualitatively the same as predicted by perturbation theory. 
Finally we have shown that the renormalized Polyakov loop defined in has the correct scaling behavior near the 
critical temperature. 

As an outlook we note that we have studied the free energy of static charges in fundamental representation. It will 
be interesting to investigate the free energy of static charges in other representations. Some work in this direction 
was done in ISll. 



Colour Singlet Correlators 


P 


T/T c 




2.3019 


1.01 


0.884(80) 


2.3077 


1.03 


1.356(52) 


2.3134 


1.05 


1.75(10) 


2.3272 


1.10 


2.23(4) 


2.3533 


1.20 


2.35(24) 


2.3776 


1.30 


2.50(6) 


2.4215 


1.50 


2.62(6) 


2.5118 


2.00 


2.37(9) 


2.8800 


6.062 


1.912(44) 


3.0230 


9.143 


1.812(36) 



TABLE II: Screening masses extracted from color singlet free energy. 



Colour Averaged Correlators 





T/T c 


lXavg{T) 

Point-C, d=l 


IT, extracted : 
Point-C, d=2 


rom 

Plane-C 


2.3019 


1.01 


0.424(20) 


< 


0.468(28) 


2.3134 


1.05 


1.000(28) 


0.544(20) 


1.024(16) 


2.3533 


1.20 


1.95(13) 


1.19(22) 


1.93(6) 


2.3776 


1.30 


2.31(16) 


1.332(88) 


2.296(64) 


2.5112 


2.00 


2.69(15) 


2.06(12) 


2.89(18) 


2.8800 


6.062 


3.03(10) 


2.32(10) 




3.0230 


9.143 


3.04(20) 


2.04(56) 





TABLE III: Screening masses extracted from color averaged free energy and from plane-plane correlators of Polyakov loops 
(see text). 
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